High-order Harmonic Generation and Dynamic Localization in a driven two-level 
system, a non-perturbative solution using the Floquet-Green formalism 



D.F. Martinez 

Max-Planck-Institut fiir Physik komplexer Systeme, Nothnitzer Strafie 38, Dresden 01187 

(^-^ , We apply the Floquet-Green operator formalism to the case of a harmonically-driven two-level 

' system. We derive exact expressions for the quasi-energies and the components of the Floquet 

' eigenstates with the use of continued fractions. We study the avoided crossings structure of the 

, quasi-energies as a function of the strength of the driving field and give an interpretation in terms 

of resonant multi-photon processes. From the Floquet eigenstates we obtain the time-evolution 
operator. Using this operator we study Dynamic Localization and High-order Harmonic Generation 
in the non-perturbative regime. 
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I. INTRODUCTION 



The driven two-level system is a particularly important model in physics, since it has proved to be very useful in 
describing many aspects of the interaction of matter with an electromagnetic field. In the field of atomic physics, 
this system was traditionally studied using the Rotating Wave Approximation (RWA)^, in which the counter- rotating 
term of the harmonic driving is neglected. This approximation is valid only for amplitudes of the driving field that 
are small compared to the energy difference between the atomic states most affected by the driving field. This was 
typical for experiments until the 1960's. Presently, technical improvements in laser technology allow experimenters 
to produce very strong electromagnetic fields (laser pulses) which can produce interesting non-linear effects in their 
interaction with matter. 

' The excitation of a spin 1/2 atom or electron in a magnetic field interacting with a microwave field is another 
^ , example where the strength of the driving field (RF radiation) can not be assumed to be small. 

Q ■ Also, due to the progress in lithographic techniques and tunable lasers in the terahertz domain, it has become 
O , possible to construct semiconductor quantum well structures driven by a strong electromagnetic field. In these 
systems, several interesting results have been obtained, including, coherent suppression of tunneling [also called 
,—1 ] Dynamic Localization (DL)j^, collapse of mini-bands in super-lattices'^, absolute negative conductance^, photon- 
J> . assisted tunneling^ and AC stark effect^. 

CO ' The driven two-level system has also been used to describe an atom moving through a Fabry-Perot cavityi^ In this 
system, the strong coupling regime is already accessible both in the microwave and in the visible regime. 

In all of the systems mentioned before, where strong fields are involved, the "counter rotating" term in the time- 
dependent part of the Hamiltonian does contribute significantly. Therefore, different theoretical approaches beyond 
■ RWA are required. One such method was first introduced by Autler and Townes^. Using Floquet 's theorem and 
\ continued fractions, they derived a general solution in order to investigate the effect of a radio frequency field on the 
"Y^ • ^-type microwave absorption line of the doublet J = 2 ^ 1 of the molecules of gaseous OCS. 

d Later, Shirley^- employed Floquet theory to reduce the solution of a periodic time-dependent Hamiltonian to the 

^ I problem of diagonalization of a time-independent matrix. This matrix was then used to calculate eigenvalues and 
I ■ transition probabilities for a driven two-level system and also to derive higher order (beyond linear RWA) analytic 
approximations valid in the case of weak driving. Later approaches to this problem, usually within the context of a 
spin system in a magnetic and a RF field, have involved the use of continued fractionsiiiiSii^, although the role played 
^ ^ by Floquet theory in those approaches (and in particular the role of the quasi-energy) is not clear, 
r* • A different kind of analytical expression for the quasi-energy of the system, which does not involve continued 
. !^ \ fractions was reported by Zhaoi^. It makes use of group theoretical arguments for its derivation and involves the 
. evaluation of complicated infinite sums. Apart from its intrinsic theoretical value, this solution does not provide a 
^ ■ practical way to study this system. 

There have been a variety of approaches used to obtain approximate analytic expressions. Most of them have been 
obtained in the high frequency limit^5ii§ii2,, for w ^ or in the strong field limit^-'^^, for v 3> hu. Here uo refers to 
the frequency of the driving, is the energy difference between the two levels of the unperturbed Hamiltonian and 
V is the strength of the driving. One of the most interesting features of a system described by these perturbative 
solutions is that at some specific values of the driving field strength (zeros of the zeroth-order Bessel function) Dynamic 
Localization (DL) occurs, which in the context of quantum wells means that despite a nonzero tunneling probability 
between the wells, a system initially prepared in one of the wells will remain in that well indefinitely. A manifestation 
of DL [also called Coherent Destruction of Tunneling (CDT)] was first found in a tight binding potential with an 
applied AC electric fieldeSS DL in quantum wells was first studied by Grofimann el ali^ 
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In a paper by Ivanov^- in which High-order Harmonic Generation (HHG) by diatomic molecular ions is considered, 
and then more recently in works by Santana et al.^^ and also Delgado and Gomez Llorente^"^, analytical expressions 
in terms of Bessel functions were found, valid in the (perturbative) regime where DL occurs. This regime can be 
characterized by the condition wo/^ ^ for any v, or ujq/uj <C \J v jfibj ^ 1. We will say that whenever any of these 
two conditions is satisfied, the system is in the Dynamic Localization Regime (DLR). When that regime is reached, 
our results converge to the known Bessel-function expressions. 

In this work we present a fully developed Floquet-Green formalism for time periodic systems, which we then use to 
derive an exact solution for the two level system with harmonic driving. This exact solution is expressed in terms of 
continued fractions, suitable to the study of the regime where neither the RWA nor any other perturbative approach 
works. Our aim is to present a whole picture of this system, of its eigenvalues and eigcnstates in the different regions 
of its parameter space and to obtain from the general solution a better understanding of the two main features of this 
system, namely. Dynamic Localization (DL) and High-Order Harmonic Generation (HHG). 

In section II we derive the Floquet-Green operator formalism and establish a connection between the poles of this 
operator and the components of the Floquet eigenstates. 

In section HI we apply this formalism to the specific case of the two-state system with harmonic driving. By 
making use of continued fractions, we obtain an expression for the quasi-energies. We study the structure of avoided- 
crossings in the quasi-energies and interpret them in terms of multi-photon processes. We also derive here the Floquet 
eigenstates (FES) of the system and study their dependence on the amplitude of the driving field for different values 
of the two-level energy difference. The time-evolution operator is constructed using the FES. 

In section IV, using the above mentioned operator, we study Dynamic Localization (DL), and in section V we study 
High-order Harmonic Generation (HHG). 

In Section VI we summarize our findings and give some concluding remarks. 



II. FLOQUET-GREEN OPERATOR FOR HARMONIC DRIVING 

The pioneering work of ShirlejiiS and Sambe^^ laid down the theoretical foundations for a complete treatment of 
time-periodic potentials, based on the same mathematical tools already developed for time-independent potentials. 
Of great importance among these tools is the Green's function, whose definition and application for time-periodic 
systems has not been clear until recently. A Floquet-Green function method for the solution of radiative electron 
scattering in a strong laser field was introduced by FaisalS^. In this work we present a fully developed picture of the 
Floquet-Green's operator formalism for general time-harmonic Hamiltonians that was introduced by the author in a 
previous work^^, and apply it to the case of a two-state system with harmonic driving. 

We start by considering a Hamiltonian of the general form: 



E(t) =Ho + 2Vcos{ujt), (1) 

were Hq and V are Hermitian operators in the Hilbert space (7i) of the system. Because of the periodicity of the 
Hamiltonian, according to Floquet's theorem, the solutions to Schrodinger's equation ih-g^\'^(t)) = H{x,t)\'i'(t)) are 
of the form 

|5-^(t)) =e-*^*/''|(/.^(t)), (2) 

where \^%t)) = \r{t+^)). 

Inserting this into Schrodinger's equation one arrives at the eigenvalue equation 

H^m'^it))^e\^'^it)), (3) 

where H^{t) is defined as 

H^{t)^Hit)^zh-l-^. (4) 

As pointed out by SambeSit, since Eq.® is an eigenvalue equation, it can be solved using the standard techniques 
developed for time-independent Hamiltonians, provided we extend the Hilbert space to include the space of time- 
periodic functions. In this extended space, the time parameter can be treated as another degree of freedom of the 
system. A similar concept is used in classical mechanics, and gives rise to the concept of a "half degree of freedom 
when dealing with time dependent Hamiltonians. 
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A suitable basis for this extended Hilbert space (7?.) is {\a) ^ |n),...}, where is a basis for the Hilbert 

space H of the system, and we define {t\n) = ^-inuit^ ^j^.]^ ^ integer. Clearly spans the vector space (T 

of periodic functions, and therefore, TZ = H^T. In this basis, Eq.Q becomes a matrix eigenvalue equation of 
infinite dimension with an infinite number of eigenvalues. It is not difhcult to proof that if is an eigenvalue, with 
corresponding eigenvector \(f)'^' {t)), then + m is also an eigenvalue (all quantities are assumed to be in units of hu), 
with corresponding eigenvector e™'^*|(/)'^' (t)). Accordingly, the eigenstate corresponding to eigenvalue 

ei + m has the same structure than the eigenstate corresponding to e^, except that it is displaced by m on the energy 
axis. Because of this, to find all the eigenvectors and eigenvalues of the Floquet Hamiltonian one needs only to consider 
— i < e < ^. We will use the letter e to refer to the Floquet eigenvalues restricted to this interval and call them 
" quasi-energies" . Clearly, any Floquet eigenvalue can be written as = Si + p for some —^<Si<^ and some 
integer p. It can be shown that, in general, there are N distinct quasi-energies (except for accidental degeneracies) if 
the Hilbert space H is N dimensional. 

This periodic structure in the eigenvalues does not mean that the " replica" eigenstates have no relevance; they are 
also valid solutions of Eq. Q and are essential for completeness in the extended Hilbert space 7^^. They also allow 
us to understand some features of the quasi-energies of the system in terms of avoided crossings between "replica" 
eigenstates. 

The Floquet-Green operator for the Floquet Hamiltonian in Eq. is defined by the equationS^, 



[iE-H^{t')]giE,t',r)^iSrit' ~t") , (5) 

where Sr{t) is the r-periodic delta function (t = ^). 

In terms of the complete (infinite) set {\<p'^' (t))} of eigenfunctions of the Floquet-Hamiltonian (Eq^J, the solution 
for Eq.® is 



I 

From the previous discussion about the eigenvalues and eigenfunctions of the Floquet Hamiltonian, we can write 
the Floquet-Green operator entirely in terms of the eigenfunctions corresponding to values between — ^ and i. We 
will denote these eigenvalues as e^. Using this, the Floquet-Green operator can be written as 

g(i^,^-,t-)^VVe-^"(*'-*")'^"^^^)^<^"^^")' , (7) 

E — — p 

where 7 = 1, ...N for TL being N-dimensional, and p = — oo, oo. 

Operating on both sides of this equation with gwnwt ^-muit ^if'^if .^^g obtain 



7,P 



where 



and 



''' Jo Jo 



e 



imujt' I 16 (j-f 



W{t'))dt'. (9) 



For m — n — Q, 



go,o(i^)=E j;,,\_p l^P^)('^P^I- (10) 



7.P 



This last equation shows us the relationship between the residue of the operator Qo,o{E) at the pole e-^ + p and 
the Fourier component of the Floquet eigenstate |<^^"'(t)). Notice that |(/)^^ (i)) = X]p ^~'^"*I'^p^)- From this we 
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conclude that, corresponding to each pole at E = e-y + p of the Floquet-Green operator Gofl{E) there is an oscillating 
term of the form e"'?"^* in the Floquet eigenstate specified by the quasi- energy e^. 

Let us go back to the task of finding an exphcit solution to the Floquet-Green operator in Eq.©. Applying 
r /o'/o'^""'^* e~™'^* '(it"(it' on both sides of this equation, and using Eq.lQ together with the definitions above we 
get 

[1[E + m) - i/o]5m,n - V[Q,n+l,n + ^m-l,n) = l<5„i,„ (11) 

The explicit solution of this equation, in terms of matrix continued fractions follows from*^. The resulting expression 
for the operator Go,q{E) is 

GQ^oiE) = [IE - Ho - V,s{E))-\ (12) 

where 

VME)^V+{E) + V-g{E) , 

with 

V,%{E)^V i J V , (13) 

E±l-Ho~ V J V 

E±2- Hq-V V 

We now specialize this result for the case of a two state system driven by a classical single frequency potential. 



III. HARMONICALLY DRIVEN TWO-LEVEL SYSTEM 



In this section we apply the formalism previously presented to the case of a driven two-level system of the form 

H = y cr, + 2vcos{ujt)a^ , (14) 

where i?, (Dq and v are given in units of fko. In the context of an electron interacting with an electric field, the amplitude 
2v corresponds to the (static) dipole moment of the electron times the applied electric field, i.e. 2v = fiEo/hcj, (the 
factor of two has been introduced for convenience) and the classical electric field is of the form E{t) = EQCOs{u!t). 
In a double quantum-dot realization of this Hamiltonian, the non-driven part, ^CTx, describes the tunneling between 
the states localized on each dot, {|1),|2)}. This tunneling gives rise to a splitting of ?iuJo = luoKuj between the 
energy levels —■^,4-^^, whose corresponding eigenstates we denote by {|a),|6)}, and where \b) = :72(l"'^^ ^ 1^^^' 
and \a) = :^(|1) — |2)). The Hamiltonian in Ea. H14|) is written in the basis {|a), \b)}. For an atomic system, ^cr^ 
describes two bound states of the atomic (or molecular) potential, or, in the context of ionizing systems it could also 
be thought of as describing a bound state and the continuum22i2S. An electron moving in a Fabry-Perot cavity^ or a 
1/2 spin system with an applied magnetic and microwave fieldii^iSiii, are also examples of systems where this model 
has been used. 

Continued fractions methods have been used before to study this system, specially in the context of spin 1/2 systems 
in a magnetic and RF field, starting with the seminal work of Autler and Townes^. Also, a previous treatment of this 
problem where the time dependent Green's function and continued fractions were used, although without reference 
to the Floquet formalism can be found in the work by Gush and Gush^°. There have been several perturbative 
approaches to this problem, as it was mentioned before. Here, we will focus on the properties of the exact solution, 
and in particular we will study its dependence on the strength of the driving field and the connection with HHG and 
DL. 

To find the exact solution to this Hamiltonian we need to evaluate Eqs. (|12I13|I . for the case V = vax and Hq — ^(Jz- 
It can be checked easily that, given the off-diagonal form of V and the diagonal form of Hq, the matrices V^{E) are 
diagonal. Its diagonal components will be denoted respectively by {v^ {E) , (E)} . The dynamic effective-potential 
Vcs{E) is therefore also diagonal and of the form 



V,siE) 



v+{E)+v-{E) 

v+{E)+v^{E) 



(15) 
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The 

^a/fc functions are given by the following coupled recursive relations 



E±l - ^ -v^{E±l) 



(16) 



Replacing into (and vice versa) we obtain the de-coupled recursive relations 



E±l 



2 E±2+^-vi{E±2) 



E±l 



g'^iE) 

g\E) 



From Ea. (|12|l . the solution for Ga.o{E) is also a diagonal matrix of the form 

GoAE) 

with 

g'^iE) = ^ 



1 



E+^-v^iE)-v-{E) 



LOq 



(17) 



E±2-^-v+{E±2) 



(18) 



g'iE) = 



E-^-v+{E)~v^{E) 



(19) 



It can easily be seen from Ea. (|17|l that this system possesses the symmetry v^{—E) — ~v^{E), which from Ea. H19|l 
gives us ga{—E) = ~gb{E). This means that to find the two quasi-energies of the system one only has to look for a 
pole in either ga{E) or gb{E) in the range —0.5 < E < 0.5; the other pole is symmetrically located on the opposite 
side of the (quasi) energy axis. 

A. Poles of GofiiE) and the quasi-energies of the system 

As discussed before, the quasi-energies of the system can be obtained from the poles in ga{E) and gb{E). We 
showed that it is sufficient to study the poles in only one of these functions since for any pole in ga{E) aX E — E* 
there is a corresponding pole in gb{E) at E = —E*. This implies the existence of two quasi-energies in the system, 
£,—e. The quasi-energy e can be obtained from the solution to gb{E)^^ = in the interval —0.5 < E < 0.5. From 
Eqs. H17|) . H19|I . the poles of gb{E) satisfy the equation. 



S + 2- 



5 + '6 + Qo 



S - 1 + Cbo — 



6-2- 



S + 4-- 



d — 3 + ujQ 



6-4 



(20) 



where S = E — The quasi-energy is obtained as e = (^ + i5) mod 1. 

In Fig. 1 we show one of the quasi-energies of the system, for different values of ujq. These functions are plotted 
as a function of the driving strength v. We have included the resonant case lLiq = 1 and several values in the high 
frequency regime where ujq < 1. In the limit when ujq << max{l, y/v), the function e{v) « ^Jo(4w). This result 
was first obtained by ShirlejilS. It is precisely in this regime where it was 

shown22-?3 |;}jat DL appears as a dominant 
feature. As ljq decreases we observe a very quick convergence to the approximate result e{v) « ^Jq{Av). Already 
for ujQ = 0.25, the difference between the exact result for e{v) and the approximation cannot be observed at the scale 
of the graph. For uiq = 0.1 we plotted both e{v) and —s{v) to make the location of the zeros of Jo(4u) more visible. 
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As we can also see in Fig. 2, increasing ojo always displaces the zeros of e{v) towards smaller values. This behavior 
was also observed by Villas-Boas et alj^ in the case of a double quantum well with an intense AC electric field and 
an applied magnetic field. Notice also that, as it is easy to derive analytically, the only case when the initial slope of 
the quasi-energy is not zero is for the resonant case.^^ 

In Fig. 2 we show the quasi-energy as a function of the driving strength v, for low frequencies, where luq > 1. The 
resonant case is included for comparison. Also, the dotted line shows the functions ±0.05Jo(4y) in order to indicate 
the location of the zeros of this Bessel function. As one increases ujq, the amplitude and sharpness of the oscillations 
around e = increases, the location of the zeros shifts to the left, and the parabola that describes (except in the 
resonant case) the initial behavior of the curve opens up. For (Dq ^ 2 one can get from Eg. 1201) that the initial behavior 
of the quasi-energy is parabolic and has the form 

Siv) « (21) 

This remains a good approximation for all values of v until 6{v) « , and provided e{v) does not reach the edges of 
the Brillouin zone (-0.5,0.5) before that. If it does so, then, the quasi-energy experiences an avoided crossing, whose 
sharpness is proportional to luq and decreases with the value of v at which the anti-crossing occurs. 



B. Avoided crossings 



In Fig. 3 we can see that the structure of the quasi-energy is determined by a regular sequence of avoided crossings 
which get broader as the strength of the driving field v is increased (resonant case shown in the figure). 

Let us consider the first avoided crossings, the ones along the line marked "3 ph. res" in Fig. 3. To zero order 
in V, the eigenvalues and eigenvectors of the Floquet-Hamiltonian are the state vectors |a)|7t), |&)|"-), for any integer 
n. In this basis, the time-dependent potential can be written as V = v(Tx{a -\- a^). When the effect of one photon 
processes (in the RWA) is included, one obtains that the new set of eigenstates is {|n,±)}, for any integer n, and 
where |n, ±) = -^(|a) I"-) ± 1^) I"- ^ !))■ It's corresponding eigenvalues are = + n±v. This one photon process is 
in a sense an avoided crossing, with a width that goes to zero as u — > 0. One can see that, near w « 1, the eigenvalues 
corresponding to the eigenstates |n, ±) and |n±2, ^) would cross if only fist order processes where allowed. To second 
order this degeneracy is not lifted since {n ± 2, =F|V^^|n, ±) = 0. Only when third order processes are included, this 
degeneracy is lifted, since it is easy to check that {n ± 2, =F|V^'^|ri±) ^ 0. Third order processes are therefore involved 
in the avoided crossings near z; « 1. 

For the resonant case, after the first avoided crossing {v « 1), for example between the states |1, — ) and | — 1,+), 
and because of considerable overlap between the 3 photon resonance and the 5 photon resonance, this states do not 
completely exchange their identities, as it would be the case in a well isolated avoided crossing. An "identity exchange" 
between Floquet states can clearly be observed for values of wq and v for which there is a sharp avoided crossing (i.e. 
(Do = 3.0 in Fig. 5). In Fig. 3, for the resonant case we can see that even though the structure of the eigenstates after 
the first avoided crossing is not necessarily of the form |n)), the points where the quasi- energies of these "first-order" 
FES cross, come close to the actual location of the avoided crossings. This is an indication that those states are 
dominant in the structure of the eigenstates, a fact that can be checked by comparing the actual structure of the FES 
in Fig. 6, at some amplitude u, with the corresponding "first order" FES eigenvalues that intersect at the closest 
avoided crossing in Fig. 3. For w « 5 the peak components in the FES are only one unit of energy further out that 
the expected value from the simple picture in Fig. 3. It is an interesting feature of this system the fact that the FES 
obtained for small driving amplitudes constitute the "bare bones" structure of the FES and to a good extent they 
determine the avoided crossings structure of the quasi-energies (at least for a larger-than-expected range of driving 
field amplitudes). 

This avoided crossings structure, to a certain approximation, can be thought of as mostly due to odd-photon 
processes between FES. Even more, we will argue that at least for a range oi v (0 < u ^ 4 for the resonant case), the 
main contribution to the avoided crossings is a three photon process for the first avoided crossing, 5 photon for the 
second avoided crossing, and so on (in the resonant case). For any value of (Do we find that the number of photons 
involved in the first resonance depends on the initial distance between the two unperturbed levels ((Do). If for example 
(Do = 3 the first avoided crossing will be the 5 photon resonance (Fig. 5); for (Do = 5 it will be the 7 photon resonance, 
etc. In general we can write that, for (Do > 1 the order of the first resonance is 2 Int( "''~"'^ ) -I- 3. The avoided crossings 
picture holds for any (Do provided v is not too big; an estimate is w ^ 4(Dq. Above this, the overlap between different 
photon processes is too big for our simple picture to apply. 

In Fig. 3. it can be seen that the states that intervene at each avoided crossing are of the form \2n + 1, — ) and 
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\2n + 1 — 2m, +), for m = 1,2, 3.... It is easy to check that 

(2n + l,-|y«|2n+ 1 - 2m,+) 7^ 0, only for q = 2m - 1,2m + 1, ... 

This is a good indication that these avoided crossings can be interpreted as muhi-photon resonances that involve 
odd numbers of photons. From Fig. 3, we see that at the avoided crossing A, from the interacting states |1, — ) and 
I — 1,+), there should be a three photon exchange between |a)|l) and |5)| — 2) and a 1 photon exchange between 
\a) \ — 1) and |6)|0). From the exact result for the components of one of the FES shown in Fig. 6 we see that the 
dominant component before the crossing is |a)|l) and therefore, at u = 1 (which is near the crossing), the two dominant 
components should be |a)|l) and \b)\ — 2); this corresponds to a three photon process. Near the avoided crossing B 
in Fig. 3, where v ^ 2, we can see from Fig. 5 that the dominant components are \a, 3) and \b, —2) which corresponds 
to a 5 photon exchange. For values of ujq and v for which the avoided crossings are sharper one gets an even clearer 
picture of those avoided crossing as resulting from the exchange of a specific odd number of photons between Floquet 
states. 

It is interesting to see how this description of avoided crossings and photon processes is reflected by the structure 
of the continued fractions in Ea. (|2()|l and how the different terms in it contribute to the features of the quasi-energy 
as a function of the driving amplitude. From examination of Eg. 120(1 it is not difficult to see that for small v, each 
term of the form S + a gives a solution near S — —a. To find accurately the solution of S for a particular value of v 
one does not need to (nor can) include an infinite number of terms in the C.F. An important question is therefore, 
how many terms to include and why? To study this, we first propose a way to truncate Ea. (|20|l so that we obtain, in 
each case, a symmetrical set of solutions. For illustration we will continue focusing on the particular case of resonance 
(ujQ = 1), as it is straight forward to apply the same ideas shown here for other values of luq. 

Ea. H20() can be put in the form 



F+{E)F-{E) ^v^, (22) 



with 



F+{E) = 6 5 , F-{E) ^ J-l + cLo 2 • (23) 

V „ V 

S + l+Qo 9 S-2 5 

+ 2 5" — 3 + 0)0 o- 



(5 + 3 + a)o 5-4 

To understand the effect of the different terms in the C.F. we will truncate the functions F^ at different levels and 
plot the solutions for S. F^ refers to the truncated version of F"^ where appears n timesS, i.e. 

F+ = 5, F+ = S-—-^ , F+ = S 



S+l+ua 



S + 2 



and likewise for F^ . 

In the different panels in Fig. 4 we plot the solutions for S{v) as obtained from the equation F^F^ = for 
ri = 0, 1,2,3 [ panels a) through d)], and only for the resonant case = 1- We also plot with a dashed line and 
for comparison, the exact result for 5 (only the two smallest roots are shown). Also, in front of each plot there is an 
insert with a diagram with lines and arrows. The lines represent each term that has been included in the continued 
fractions F^ and F~ . As we mentioned before, each term of the form 5 + a contributes with an eigenvalue curve that 
originates aX 5 — —a. There are two kinds of terms of that form in the C.F.'s of Ea. ((2()ll : 5±2n and (5 + cJo ± (2m + 1). 
The first ones contribute, at u = 0, with a solution at 5 = +2n and therefore a pole aX E = E^^ 2n. The second kind 
gives a solution for 5 = — (Dq T (2m + 1) which contributes with a pole at i? = i?6 — T (2m + 1) = £'a + (2m + 1). 
To the lines on the left we associate the states |a, 2n + 1) . We will label this lines as a + 2n + 1, for some n; the lines 
on the right we associate with the particle-field states \b, 2n) and label them as 6 + 2n. A similar picture would result 
from studying the poles of ga{E). 

When two such terms in the CF are included, then we see, in the corresponding panel in Fig. 4, a pair of solutions, 
which at u = converge to 5 = 2n and split for u > 0. We know, from previous discussions, that the upper one should 
correspond to the FES, \2n + 1, +) = \a, 2n + 1) + \h, 2n), and the lower one to \2n + 1, -) = \a, 2n + 1) - \h, 2n). 
The arrows at the end of the lines show the direction in which the corresponding eigenvalue would move when v is 
increased. The repulsion between the eigenstates |2n + 1,±) is due to the exchange of one photon between \a) and 
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1 6). This is the main process at work in panels a) and b). In panel c) two new processes enter the picture: As the 
levels move apart due to the one photon resonance, another resonance (A) can now occur when the levels |3, — ) and 
|1, +) mostly due to the transfer of three photons between |a, 3) and \b,0). This process appears for the first time for 
when the term |a, 3) is included in the CF. A replica of this three photon resonance occurs also between levels 
and 1-1,+). 

When \a, 3) and \b, —2) is included, the second process that becomes available is the five photon resonance between 
levels |3, — ) and | — 1, +) (after this level have interacted with |1, +) and |1, — ) respectively). This resonance becomes 
more clearly defined in panel d) where more levels are included. Even though, as compared to c), no new resonances 
appear in d), the inclusion of levels |o, —3) and \b, 4) contributes to sharpen the resonances already visible in c) and 
in particular the 5 photon resonance B. This shows that, as mentioned before, there are several odd-photon processes 
that contribute to a particular avoided crossing. There is however, an odd number of photons that seems to contribute 
the most and without which, if the corresponding levels are not included, the corresponding avoided crossing does 
not seem to occur. That specific number of photons is shown on the top of Fig. 3. For v > 4 this pictures begins to 
falter, and the avoided crossings begin to depend more on a higher number of photons than what would be predicted 
from Fig. 3. For v = 5 for example, a better description of the crossing would be for a 15 photon resonant process, 
not the expected 13. For even higher values of ujq this picture holds on a bigger range of the driving field amplitude, 
namely, for v ^ 4cZ'g. As mentioned before, one can see that in Fig. 5, for ujq — 3, the first avoided crossing is clearly 
a 5 photon resonance. 

As an empirical rule, which agrees with the above observations, we found that to obtain an accurate result for the 
quasi-energy up to a value v of the driving potential one needs to use the function F^^{E), with n > Int(2w). 



C. Components of the Floquet EigenStates 

In the previous subsection we showed how to obtain the diagonal part of the Floquet-Green operator and discussed 
its pole structure and their avoided crossings. In this part we study the Floquet eigenstates and the general behavior 
of them as the amplitude of the driving field changes. In Appendix A we show that the FES of this system can be 
written in the form 

\r%t)) ^ K{t)\a) + Z{t)\b) 

and 

|rW)=-Z(t)*|a) + Kitr\b), (24) 

with 

oo 

K(t)=N[ e-2»-*i^2„(-e)] , 

(25) 

n— — oo 

TV is a normalization constant and |iir(t)p + |Z(i:)p — 1. In appendix A we also show the explicit form of the 
components Z2n+ij in Ec ljCll) . It can there be seen that all components and Z2n-i are real, and therefore 
K{0) and Z{0) are also real numbers. 

In Fig. 6 we show the components of one of the Floquet eigenstates of the system, [^"^(i)), for the resonant case, 
and for eight different values of the driving field amplitude, In the upper figure of each panel we show the quantities 
Zn,Kn (with n odd for Zn and n even for Kn), and in the lower panels we show the quantities ln|^„|, Inji^nj. 

A noticeable qualitative difference between the behavior of Z2n-i and K2n is the fact that for stronger values of 
V, the former becomes more antisymmetric while the last one becomes more symmetrical around n = 0. We will 
comment shortly about the consequence that this has for the dynamic localization of the system. The fact that the 
eigenstates of the system have components that are harmonics of the driving field frequency w, translates, as we will 
study in the next section, into a (coherent) emission spectrum of the system which contains also many frequencies 
that are multiples of the frequency of the driving field. The possibility of exciting a system at a particular frequency 
with a strong field, and producing radiation that includes many harmonics of the driving frequency (in some cases 
even hundreds of them) is called High Harmonic Generation (HHG) and has been studied in many experiments^ 
Driven two-level systems have been studied in the past to understand the general features of HHG such as the 
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existence of a plateau in the spectrum followed by a sharp cut-off where the amplitudes of the harmonics decrease 

^^p^j^.. 16.21.28.35.36.37.38 

As one can see directly from Fig. 6 (resonant case), in most of the cases shown, the Floquet eigenstates show 
a particular structure with two clearly distinguishable regions. The first region which, in a non-rigorous way we 
call the "chaotic" one, can be described as corresponding to the Floquet components with \n\ ^ — 2v ~ 1 (for 
V > 1 in the resonant case). The amplitude of these Floquet components depends strongly on n and on v, and 
there are frequent sign changes. This gives rise to the "plateau" region typical in HHG. The second region, or the 
"regular" one corresponds to the components with |n| > N^, and it is characterized by having amplitudes that decay 
exponentially with n. Also, notice that the components Kn and Z„ do not change sign in that region, with > 0, 
and sg(Z2„_i) = sg(n) in the " regular" region (sg(x) = a;/|a;|). For other values of Cjq, a plateau region was always 
found whenever v ^ y/ujQ (for ujq > 1), in agreement with a similar expression obtained by KaplanSS in the limit 
ujQ >> 1. For Qq < 1, the plateau forms for any v ^ 1. 

The expression — 2v — 1 can be derived in a simple way. First, we look at the conditions to obtain, for n > 
a perfectly flat, triple plateau in the components of the FES. This means that we will assume Kn+2 — Zn+i = Kn 
(which gives three consecutive components with the same amplitude), for some n even. This condition can be achieved 
if, from Eq.jnU, {e + n) = v'^ {e + n + 1) = v. From Ea. lfTC|l we get the condition, 

n = 2v — 1 + — — e. 

If instead, we had chosen the triple plateau to be of the form Zn+2 — Kn+i = Zn ( for n odd), then, the equation for 
n would be 

n = 2v — 1 e. 

2 

If we now apply the same reasoning for a triple cusp to occur in the negative energy components, beginning at the 
component — n, we get the conditions (for n odd and even) 

71 = 2i; - 1 + y + e. 



n = 2v — 1 — + e. 

Clearly, since in general this 4 equations can not be (nor we care for them to be) satisfied simultaneously, we look for a 
condition that will produce a simple (possibly) double cusp, both at Ne and —Ng. This would describe the location of 
the two broad peaks that are, in all cases, at the limit between the "chaotic" and the "regular" region. For obtaining 
this symmetric flat cusp condition we pick the equation that is right in the middle of the set of 4 equations given 
above. This equation is clearly 

Ne = 2v-1, (26) 

with the nice property that the e and the luq dependence drops out. This is a result that should work well whenever 
there is a well defined "plateau" region, that is, for v > a/Sq when luq > 1, or v > 1 when cjo < 1- As can be seen 
in Fig. 6, for (Dq = 1 and v > I, this equation seems to describe well the location of the peaks in the spectrum of 
an eigenstate that separate the plateau from the exponential decay region. For other values of a)o it gives also good 
agreement with the data. 

To investigate further the structure of the FES of the system, we use Ea. H24|l . and evaluate a.t t — 0, 

ir^(O)) = cose\a) + sin0|6), 

|(/) ^(0)) = -sin 6*10) + cos 6*16), ^ ' 

and where, using K{0)^ + Z{0)^ = 1, we have defined the angle 9 so that K{0) — cos 6, and Z{0) ~ sin 9. 

We now concentrate on the study of the quantities cos 9{v), sin d{v), d{v), which, as we will see in the next 
sections, play an important part in the evolution of the system, and correspondingly in the Dynamical Localization 
and Harmonic Generation of the system. In what follows we will study them for three different values of the parameter 
ojq: Driving frequency above resonance, for ujq = 0.5; resonant case for — 1 and driving frequency below resonance, 
for (Do — 2. 

In Figs. 7,8,9, we show the results for wq — 0.5, luq — 1.0 and wq = 2.0 respectively. The upper panels show 
the dependence of the functions cos'^{29{v)) and sin^(26'(u)), and the lower panel the functions 9{v) and e{v). As a 
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function of the driving field amplitude v, these functions show an oscillatory behavior, which can be better appreciated 
in the lower panels, where the function 9{v) is plotted along side with s(v). An interesting feature of this functions is 
the fact that in general, and to a good approximation, there seems to be a phase shift of w tt/2 between e{v) and 6{v), 
which corresponds to a distance of Av « 0.4 between their maxima. This means that to (almost) every maxima of 
9(v) corresponds a zero of e{v) (the DL points) and vice versa, with the only exceptions occurring at the first maxima 
or the first zero of this functions. Notice also that d{v) does not oscillate symmetrically around zero, with the clear 
consequence that the peaks in the functions cos^(20(w)) and sin^(20(w)) have uneven heights, with alternating high 
and low peaks, although with an overall decrease in height. This alternating pattern of maxima seems to be a general 
feature of the system, which we have found for all the values of (Do that we examined. 

The DL regime, which can be characterized by the parameter a' = Uq min(l, l/y^) << 1, has been studied in 
several works. In this regime it was foundiSiSS that, at f = 0, the eigenstates of the system are of the form 

\r'm = |a) + CjojHoiiv)\b), 

n (28) 
10^(0)) =-[oo^Ho{Ma) + \b) 

where Hq{4:v) is the Struve function?i defined as 

n=0 

In Figs. 7,8,9, in the lower panel we have plotted with a dash-dot line the function QojHq{4:v). As it was expected, 
since in the DL regime sin 9 « 1, then sin 9 9, and 9{v) — > LUojHo{iv). The convergence of 9{v) to the Struve 
function is clearly slower for the bigger values of ljq. It is however rather remarkable that the Struve function seems 
to approximate well the function 9{v) much before the DL regime condition is satisfied {a <^ 1); for example, for 
luq = 0.5, and v = 1, or luq = 1 and w = 4, which both give a = 0.5 ^ 1 we already see a good agreement between 
those two functions. One can also see from Fig. 7, 8, 9 that the Struve function correctly gives the height of the maxima 
in 9(v), with the difference between those two functions being mostly a phase difference that tends to zero as v oo. 
From the asymptotic form of the Struve function in the limit u — s- cx) we obtain that, in that limit. 



TT / 1 

9{v) ^ Cjo — \ - — sin(4w — 7r/4) « a sin(4w — 7r/4) for v ^ oo. (30) 
4 V 27ru 



where we defined a = y^r/M a' = for v > 1. We will call a the "DL parameter". Eq. H30|l establishes a 

nice connection between the DL parameter and the amplitude of the oscillations in the function 9{v). 

D. Time-evolution operator 

To find out the time evolution of an initial state of the form \a) or \b), we can use the eigenstates at time t=0, as 
given by Ea. (|27|l . From that equation we can express \a) and \b) in the form 

\a) = cos 9\r%0)) - sin0|0^(O)), 
\b) =sm9\(j)-'{0)) + cos 9\(l)''{0)). 

From these, it is now straight forward to construct two orthogonal solutions to Schrodinger's equation which, at i = 0, 
correspond to each one of the eigenstates of Hq : 

\a{t)) = e*^*cos 9\<j)-''{t)) - e-*^*sin 6l|0^(i)) 

and 

\b{t)) = e*^*sin 9\(f>-~'{t)) + e-^^*cos 9\(p''{t)), (32) 
where |a(0)) = \a) and \b{0)) = \b). Using Ea. l24|l into above expressions we get 

\a{t)) = [cos 9K{t)e''' + sin 9Z* {t)e-''']\a) + [cos 9Z{t)e''' - sin 9K* {t)e-''*]\b) 
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and 



|6(t)) ^ [sin eK{t)e'^* ~ cos 9Z* {t)e-''*]\a) + [sin 9Z{t)e"'* + cos 9K*it)e- 



(33) 



The time-evolution operator for this system, which is defined by the equation |^'(t)) — f/(<, 0)||5'(0)), can be 
extracted, in the basis {\a), \ b)}, directly from the expressions above 



Uabit,0) 



Uaa{t) ~ul^{t) 
Uba{t) Ul^{t) 



(34) 



where, 



Uaa{t) = cos 0K{t)e'''* + sin 0Z*{t)e- 



Uba{t) = cos eZ{t)e'^^ - sin 6'i^*(t)e"*^*. (35) 
The subindex ab in Uab is there to remind us of the basis used to represent the time-evolution operator. 



Another way to arrive at the same result is by making use of the general expression [Eq.(27) in Milfeld and Wyat 



^1 



[/(t,0) = $(i)e"'^*$~^(0) (36) 

where the columns of the Floquet matrix $(t) are the components of the FES of the system, and A is a diagonal 
matrix whose components are the quasi-energies of the system. For our system, from Eq. (|24() , the Floquet matrix is 



m 
m 



~z 

K 



it) 



(37) 



and A = — ecTr 



IV. DYNAMIC LOCALIZATION 



It is of particular interest to study the situation where the the system's initial state is a superposition of the 
eigenstates of Hq. In particular, specially in the context of quantum wells, a particularly interesting initial state is one 
in which the electron is localized in one of the wells. Such an initial state corresponds to [^'(0)) ~ |1) = "^(k) + \b)), 

or [^'(0)) = |2) = :^(|a) — \b)). In this case, we are interested in calculating the probability to find the system at 
that same location after a time t. The solution can be found by applying the transformation V = -^{<^z + <^x), which 
corresponds to the change of basis {|a), \b)} — > {|1), |2)}, to the evolution matrix in Ea. (|34|l . The resulting matrix, 
Ui2 = VUabV^^, can be written as 



C/i2(t, 0) 



Unit) -u*2i{t) 
U2lit) Mn(t) 



(38) 



with 



Wii = Reiuaa) + ilmiuba), 



U21 = Re{uba) - ilm{uaa)- (39) 

It is a well-known fact that dynamic localization occurs only at the points where the quasi-energies vanish, i.e. 
e{v) = 0. At this particular values of driving field amplitude v, as it can be seen from Eas. ljSUl and ifSSll . the dynamics 
of the system becomes strictly periodic, with period 27r/cj, for any initial condition. The points where sin 6{v) — 
also produce a periodic evolution of the system (except for an overall time-dependent phase factor), but only for the 
particular initial conditions |^'(0)) = |a) or |4'(0)) = \b). 

To study the system at the DL points, we will look at the quantity |uii(i)|^, which gives the probability (as a 
function of time), that the system, if localized in a well (e.g. |1) ) at time t — Q, would remain in that well at a later 
time t. For e{v) = 0, and from Eas. (|35l39|l we get 



ux^i{t) = cos e[Re{K{t)) + iIm{Z{t))] + sin e[Re{Z{t)) + iIm{K{t))] , at DL points. 



(40) 
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Since this probability is periodic, mi.i(0) = ui^i{2'k/uj) = 1, it is interesting to see what is the value of this function 
half-way through a cycle of the field, that is, for utt = tt. From Ea. H6U|l we can deduce that K{'k/uj) = K{0) = cos 9, 
and Z(tx jiS) = — ^(0) = —sin Q. Using this, we get from the above equation 

u\^\{jxluj) — COS 6^ — sin 8^ — cos(20) , at DL points. (41) 

We will show in the following figures that the function cos(26') plays an important role in quantifying how good is 
the localization of the system at the DL points. 

In the lower panels of Figs. 7, 8, 9, we have placed dots on the quasi-energy functions e{v), at different values of v. 
In the upper panel in Figs. 10,11,12 we show, for each one of those points, the corresponding probability as 
given by Eqs. H35|) . H39|) . This probability gives the location of the particle, as a function of time, when starting from 
the initial state |1). In Fig. 10 we study the case ujq = 0.5. In this plot one can see that the probability evolves in a 
manner that is typical of the high frequency regime which has been amply studied in the literature. The function can 
be described as of the form cos{et)+ small time-periodic oscillations. Clearly, at the DL points (A,B,C,D), when e = 0, 
the probability is exactly periodic (period 2ttll!) and oscillates between 1 and a value that depends on the driving field 
amplitude v and on a)o- As derived before, in the middle of the driving period, the probability = cos^{29{v)). 

Given the correspondence between the maxima of 9{v) and the zeros of s{v), we can conclude that the minima of the 
function cos^{26{v)) tells us the value of the probability |ui_i(t)p half way through the driving field period {t = tt/u;), 
at the corresponding DL point (with the exception of point A). This value of |Mi_i(i)p is very important because it 
gives us information about the extend to which the particle is localized when the DL points are reached. As we can see 
in the lower panels of Figs. 10,11,12, when the function |Mi,i(t)p has a local minimum at t — tt/lu, then cos^(26'(u)) 
is close to be the absolute minimum value that |Mi_i(t)p reaches through out a whole period, Pi = r7im(|ui^i(t)p). 
It is therefore an important number for describing the amount of localization that can be achieved in the system. 
When the fimction |7Xi_i(i)p has a local maximum at t — tt/uj, then cos^{29{v)) is still equal to that value but this 
is not necessarily close to Pi (specially when we are far from the DL regime). This is the case at points B and D in 
those figures. Surprisingly, however. Pi, for cases B and D, is given, to a certain approximation, by the value of the 
function cos^(20(i;)) at the minima that corresponds to the next DL point. We do not know why this should be so 
but it was found to be the case in all the cases that we examined. According to this, the DL points can be grouped 
in pairs, with (B,C) having a similar amount of localization, given by cos^{29{vc))- points (D,E) with a localization 
given by cos'^{29{ve)), and so on. It is therefore the deeper minima in cos^(2^^(ti)) that give us directly the amomit 
of localization in the system. 

An interesting question to ask, because of practical applications, is whether there is a simple relationship between 
the amount of localization P;, ujq and v. In other words, given a particular amount of localization P/, and a given 
(Do, what is the required driving field amplitude for the system to achieve such level of localization? To answer 
this question in the DL regime (were we expect a high level of localization, P/ < 1), one can use Ea. (|30|l to derive 
an approximate relationship between P/, cDq and v. At the maxima of 9(v) (Dl points), 9{v) ~ a, which gives 

Pi = cos(26l) « 1 - 46*2 = 1 - l^^i, from which one finally gets 

(42) 



8(1 -P) 

This simple equation, which has been derived assuming v > 1, a << 1, already gives good results even for a 
localization as low as P/ = 0.5, with the accuracy improving further for higher values of P;. For P/ = 0.5 the result is 
V ^ 0.78 <Dg, for Pi = 0.9 one obtains v ^ 4 and for P; = 0.99 it gives d ^ 40 Wg. 

V. HIGH-ORDER HARMONIC GENERATION 

The expected value of the dipole moment of the system can be calculated from the expression 

{d{t)) = {^it)\X\^{t)) = M(*(0)|a.(<)|*(0)) , (43) 

where 

a,{t) = Uab{t)^<J^Uab{t), (44) 

and using Eas. (|5Hl . we get 



ReKait)Uba{t)] Ka{t)-uUt)r 

<ait)-ulit) -Re[uUt)uUt)]- 



(45) 
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For a general initial condition of the form \'i'{0)) = Ca\a) + Cb\b) we get 

{d{t))/^, = Re[uUt)uta{m\ca\' - IcfcH + 2i?e [c,c^ (^L (t) - uUt))]- (46) 
If we now define the following time-periodic functions, 

R{t) = Re{K{t))Re{Z{t)) + Im{K{t))Im{Z{t)) (odd harmonics), 

Q{t) = Re{K{t))Im{K{t)) - Re{Z{t))Im{Z{t)) (even harmonics), (47) 
P(t) = Im{K{t)f + Re{Z{t)f (even harmonics), 

then, using this, and plugging Eas. H35|) into Eg . (1461) . we obtain (after some lengthy algebra) the final expression 

{d{t))/^i = [-sin(26')cos(2et) + 2cos{29)R(t) + 2sin{20)sm{2et)Q{t) + 2sin(26l)cos(2ei)P(i)] (|cap - |cf,p) 

+ [cos(26l)cos(2et) + 2sin(26i)i?(t) - 2cos(26i)sin(2ei)Q(t) - 2cos(26l)cos(2et)P(t)] Re{2cacl) (48) 
+ [-sin(2et) - 2cos{2et)Q{t) + 2sm{2et)P{t)] Im{2cacl) . 

The quantities 0, e, P{t), Q{t), R(t) depend on v, even though, we do not explicitly indicate so in this notation. 
We will refer to the initial condition where |cq| = 1 or \cb\ — 1 (and therefore CqcJ = 0), as the "Optical" Initial 
Condition (OIC); the case Ca — as the "Tunneling" initial condition (TIC) and the case Ca — iicb as the "Complex 
Tunneling" Initial Condition (CTIN). 

From the expectation value of the dipole moment, one can obtain the (coherent) emissions spectrum of the system 

2 

from the expression \d{n)\^ = i J^'"^'^ dte'^^* {d{t)) . 

Ea. l|48|) deserves many comments. Notice that generally, there is low- frequency generation, with frequency 2e for 
all initial conditions. Only at some particular driving amplitudes this can be prevented. For OIC the low-frequency 
component can be suppressed only when sm{29) = 0; for TIC only when cos(20) = 0. There is always low-frequency 
generation at any driving amplitude for CTIC. There is odd- harmonics generation, at frequencies (2n — l)uj and due 
to the terms with R{t). This harmonics are generated except when cos(26') = for OIC and except when sin(20) = 
for TIC. Interestingly enough, CTIN can not generate odd harmonics. Even harmonics, at frequencies 2nuj ± 2e occur 
because of the presence of terms (sin(2ei) or cos{2et))Q{t) and (sin(2et) or cos{2et))P{t), and are generated except 
when sm29{v) — for OIC and except when cos20(w) = for TIC. They are always present for CTIC. 

For a general driven Hamiltonian that possesses the Dynamical Symmetry (DS) H{x,t) — H{—x,t + t/2) it has 
been proven^ that the HHG from a Floquet eigenstate consists of only odd-harmonics. Our Hamiltonian, Eq. 114|l . 
is certainly of this kind since the quantum operator that corresponds to x is smd therefore, if — cTj,, and 

t ^ t + t/2, the Hamiltonian in Eq. H14|) remains the same. According to this, a good check for our result can be 
done by replacing in Ea. (|48|l the initial condition Cq = cos 9, Cb = sin 9 which corresponds to the initial state being 
already a Floquet eigenstate [see Ea. H27|l ]. After doing that, one obtains {d{t)) / jj. = 2R{t), which means that this 
state only generates odd-harmonics, in agreement with the dynamical symmetry arguments developed in the citation 
above. 

As our result (|48l) shows, in a two-level system it is possible to generate both even and odd harmonics. This has 
also been found in numerical calculations for a double quantum well**^ . The kind of harmonics present in the spectrum 
is dependent on the initial condition, the driving amplitude and also the phase of the driving field. In this work we 
have chosen a driving force of the form cos(a;i); however, in the most general case of time dependence of the form 
cos(wt -|- (j)), the driving phase i/i, in general, plays an important role in determining the harmonic content of the 
spectrum. For the case sin(wt) see Gauthey et al^. 

It has been argued^ that in systems that get ionized during the course of the interaction with the field, only one 
resonant Floquet state (the one with the longest life-time) contributes to the harmonic spectrum and therefore this 
could explain the absence of even-harmonic generation in most experiments. There are several experimental results 
that involve lasers interacting with solid targets, in which even-harmonics and odd harmonics are generated. Even 
though ionization takes place, the presence of even harmonics has been attributed to the lack of DS in the Hamiltonian. 
For systems where the main source of harmonic radiation are bound-bound transitions, it is not clear how to justify 
the use of only one of the FES in the calculation of the HG of the system. In those kinds of system we will argue 
in future work that, if the evolution of the system is adiabatic under the effect of a laser pulse, then again, only 
odd-harmonics are generated. For non-adiabatic evolution, we expect to see generation of even and odd harmonics 
and also dependence on the initial phase of the driving. It is interesting to check our results for the time dependent 
dipole moment of the system in the DL regime. As we showed before, in the limit of perfect DL, where ^ cx) or 
(Do 0, one gets 9 = 0, K{t) becomes pure real, Z{t) becomes pure imaginary, and therefore R{t) = Q{t) = P{t) = 0. 
The result for the dipole moment is therefore 

{d(t)) 111 = coii{2et)Re{2caC*b) - sin(2ei)/m(2c„c^) = 2|caCb|cos(2rf -I- 7), (49) 
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with 7 being the phase difference between Ca and Cb- This is the zeroth-order expression in the first-order perturbative 
results obtained by Delgado and Gomez-Llorente^^. Clearly, in this limit there is no HG. From this it is clear that in 
the case of a driven two- level system, DL and HHG are at odds; the greater the localization, the smaller the HHG. 
The amplitude of any harmonic component of the emission spectrum of the system (except the low frequency one) , as 
a function of increasing driving field strength v should initially increase, then plateau (with oscillations) and finally 
decrease as v is further increased. 

Near the DL regime we can obtain some expressions and compare them with the perturbative results obtained in 
other worksiSiSLS^. For v > 1 and a = ^/^^, we found that the behavior of the functions K{t) and Z{t) is 

K{t)^Re{K(t))+iati{t) 

Z{t)^aC{t)+iIm{Z{t)). ^ ' 

Here K(t) and C(i) are periodic functions (even and odd respectively) with amplitude close to one. As we have shown 
before, in this regime, the amplitude of the oscillations in the function 0{v) is also proportional to a. According to 
this, and using Ea. (|47|l . we see that in this regime, R{t) cx a, Q{t) cx a, P{t) cx . We now rewrite Ea. H48|) in a form 
that makes evident the magnitude of the different terms that contribute to the dipole moment, 

{d{t)) w [-a^cos(2ei) -I- aR{t) + 2a^hin{2et)Q{t) + 2a^ecos{2et)P{t)\ (|cap - |c6p) 

cos(2et) - 2as\n{2et)Q{t) + Aa^ORit) - 2a^cos{2et)P{t) Re{2ca4) (51) 
-sin(2e<) - 2acos(2et)Q(t) + 2a^sin(2ei)F(t) Im{2cacl) . 

where R{t) = R{t)/a, Q{t) = Q{t)/a, Pit) = P{t)/a'^, and 6 = 9/a ^ sin(4w — 7r/4). From this expression we can 
see that in the DL regime, QIC generates low frequency and odd-harmonic radiation to second order in a and even 
harmonic generation is forth order in this parameter. For TIC there is strong low frequency generation independent 
of any of the parameters of the system (provided we are in DL regime ); even harmonic generation is second order and 
odd harmonic generation is forth order in a. For CTIC the situation is similar to TIC except that no odd- harmonic 
radiation is generated at all. 

As mentioned before, from the general expression for the dipole moment, since e = e{v) and 9 = d{v), with 
successive points where either 9(v) or e{v) vanish, there is a lot of variation in the composition (in terms of odd or 
even harmonics) of the spectrum of this system as a function of the amplitude of the driving field. The interesting 
exception being the CTIC where even harmonics are always generated and odd harmonics are never generated. 

In Fig. 13, 14 and 15 we show the emission spectrum for three distinct initial conditions (OIC, TIC, CTIC) and 
for two values of the driving field amplitude. What is common to all this plots is the typical profile of harmonic 
generation that has been found theoretically and experimentally in different systems. It consists of a plateau where 
the amplitudes are of similar magnitude (fluctuations are of one to two orders of magnitude typically), followed by a 
frequency cut-off after which the harmonic amplitudes decay very quickly (exponentially in ulo). 

From the structure of the FES, like the one shown in Fig. 6, one can see that there is a cut-off in the harmonic 
components of the eigenstates, which occurs at = 2v — 1 (for v > 1). The location of the cut-off in the emission 
spectrum of the system follows from the location of the cut-off in the eigenstates, since the functions R{t), Q{t), P{t) 
are second order in the FES components. R{t), Q{t), P{t) therefore have cut-offs at Ncut — 2A^e — Av — 2. This is the 
same cut-off that we expect to see in the emission spectrum of the system. 

As we can see in Fig. 13, 14, 15, this result describes quite well the location of the last peak before the components 
start to decay exponentially. The linear dependence of Ncut on the amplitude of the field implies that Ncut oc -s/T, 
where I is the intensity of the field. This result is expected for HHG in systems where bound-bound transitions are 
dominanli^Si^. For dominant bound-continuum transitions^^, Ncut oc /. 

In Fig. 13 we show the emission spectrum for OIC, and for two different values of the driving amplitude. For 
V — 4.15 (which gives 9 ~ 0), we see, from the first line of Ea. (|51|l that one expects no low frequency component, odd 
harmonics and no even harmonics. For v — 4.5 (which gives £ ~ 0), we obtain a strong low frequency component, 
strong odd- harmonic components and very weak even-harmonics components that are due to the function P{t) only 
(since at e = 0, sin(2e<:) = and therefore Q{t) does not contribute to the spectrum). This is a special case since in 
general, in this regime, Q(i)'s contribution to |d(il)p is two orders of magnitude bigger than the one from P{t). 

In Fig. 14 we show the emissions spectrum for TIC for the same values of v as before. In the 9 — case, we can 
see, from the second line in Eg. 1)5 1(1 . that there is a strong low frequency component (O = ±2e), no odd- harmonics 
and split even harmonics {il — 2nzL 2e) . For the case e = there is a zero frequency component and also odd and 
even harmonics (not split) of the same magnitude. 
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Finally, in Fig. 15 we show the case when the initial condition is CTIC. From the third line in Ea. (|51() we see that, 
for the case 9 = we should get a strong low frequency component and split even harmonics. For e = case there is 
no static component, only even harmonics (not split). 

For the characteristic plateau of HHG to appear, it was shown^^ that v > y/uJo- In this region we have found 
numerically that cos(20) ^ 0, which means that in the OIC, overall, the odd harmonics are much more stable with 
respect to changes in the driving amplitude that the even harmonics, which as a whole, do fluctuate with v (following 
the behavior of sin(20)). As one approaches the DL regime the odd harmonics will dominate the spectrum for the OIC 
(which is the initial condition most easily obtained experimentally) . It is easy to derive an estimate for the conditions 
under which odd-harmonics will dominate over even harmonics in the OIC. For that we take a > 0.1, so that there is 
at least one order of magnitude difference in their contributions to Eq.(49). From a = V^"^ one gets the estimate 

When this condition is satisfied, odd harmonics will clearly dominate the spectrum in the OIC. 



VI. SUMMARY AND CONCLUSIONS 



In this work we used a complete Floquet-Green operator formalism for the analytical solution of a harmonically 
driven two level system. From this operator we were able to completely solve the system and obtain the quasi-energies 
and the eigenstates in terms of continued fractions. 

We found that the oscillatory behavior of the quasi-energies of the system can be interpreted as due to multi- 
photon resonances that occur in an increasingly periodic way as the strength of the driving field is increased. This 
multi-photon resonances involve an odd number of photons and produce avoided crossings centered at the edges of 
the energy Brillouin zone e = ±0.5. The width and amplitude of them increases with the driving field amplitude. 
The structure of the continued fractions from which the quasi-energies can be obtained made it clear that this is the 
case, since each term in the continued fraction was found to be related to the inclusion of a particular basis state in 
the calculations, and each avoided crossing studied did only appear when the corresponding term in the continued 
fraction was included. 

We were able to obtain analytical expressions for the components of the Floquet eigenstates and from them we 
constructed the time-evolution operator of this system. The structure of the eigenstates revealed a characteristic 
pattern, with a "chaotic" region where the components are sensitive to small changes in change sign and have 
zeros, and a "regular" region where the components depend less sensitively on v, and decay exponentially with 
the frequency nuj. The cut-off frequency that divides the two regions and gives the position of the largest floquet 
component of the eigenstates is proportional to the driving field strength v, Nf, = 2ti — 1, a result that is valid for any 
Wo (provided v > VuJq)- 

From the time evolution operator we were able to study the Dynamic Localization phenomena which is one of the 
important features of this Hamiltonian. It is known that DL is obtained at values of the driving field strength for 
which the quasi-energies vanish. In Quantum Computation, this localization mechanism can be of great importance 
for controlling the state of a q-bit (as realized for example in a solid state devise such as a double quantum well or 
a Josephson Junction), since the amount of tunneling of an electron between two adjacent wells can be varied with 
the intensity of the driving field. In this work we obtained a correlation between the DL points and the maxima of 
the quantity 9{v) that describes the Floquet eigenstates of the system. Also, the value of the function cos^(26'(w)), 
at its lowest minima, was found to give a good approximation to the degree of localization of the system at the 
corresponding DL points. These points appear in pairs with simmilar degrees of localization corresponding to each 
pair. We also found the equation 



8(1 -PO' 

which, at the DL points, relates the degree of localization Pi, with the energy level difference loq, and the amplitude 
V of the driving field (all quantities in units of Hlu). This equation can be used to easily estimate the amount of 
localization achievable in this system for a given set of parameters. 

For the study of High Harmonic Generation in this system we calculated the expected value of the dipole moment 
which revealed that the emission spectrum, in general, can contain a low frequency component and odd and even 
harmonics. The particular components and their amplitudes depend strongly on the initial state of the system and 
on the specific driving field amplitude. We found that there exists an initial condition (Complex Tunneling) of the 
driven two-level system for which no odd-harmonics can be generated at any driving field strength. For each one 
of the "pure" initial conditions, namely Optical (OIC), Tunneling (TIC) or Complex Tunneling (CTIC), the driving 
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amplitudes where the 9{v) = produce either pure even harmonic generation or pure odd harmonic generation, 
depending on the particular initial condition. For all other driving amplitudes there will be in general generation of 
both kinds of harmonics (except for CT initial condition). When the initial state is a Floquet eigenstate, we found that 
only odd-harmonics are generated, a result that is expected from the Dynamical Symmetry (DS) of the Hamiltonian. 
In general we can say that the emission spectrum of the driven two-level system is sensitive to different parameters, 
such as the initial preparation, the amplitude of the driving and its phase (phase dependence will be explored in 
future work). 

For small values of a, but still not deep into the DL regime {a w 0.15) we showed the different magnitude of the 
terms that give low frequency, even harmonics and odd harmonics. In the Optical Initial Condition (OIC) we found 
that odd harmonic generation is stronger, with even harmonics weaker and even nonexistent for particular values of 
the driving field amplitude (when 6(v) = 0). For Tunneling Initial Condition (TIC), even harmonics dominate, except 
at the DL points (e = 0), where even and odd harmonics have approximately the same amplitude. For Complex 
Tunneling Initial Condition (CTIC) no odd harmonics are generated. 

Well inside the DL regime, odd harmonics are almost exclusively generated in the Optical Initial Condition, and 
even harmonics in both of the Tunneling Initial Conditions. In this regime of high localization. Harmonic generation 
decreases with decreasing a, a situation that clearly places both effects on opposite sides of the parameter space, 
although with significant overlap between them. In general, for a driven two level system, the strongest HHG should 
occur for higher values of c^o and for driving amplitudes v ^ lOwg ( for a > 0.1). For strong DL, smaller values of Qq 
and/or bigger values of v are required, v > IOQq (for a < 0.1). Strong HHG in a two-level system is an eminently 
non-linear, non-perturbative phenomenon. 

The driven two-level model has been shown to be useful in understanding the basic features of high-harmonic 
generation experimentally observed in atoms and molecules, even though obviously the full complexity of the general 
processes in atoms and molecules can not be accounted for with this model. The model predicts an emission spectrum 
with a plateau and a cut-off frequency {N^ut = 4i; — 2), whose location depends on the square-root of the intensity 
of the driving field (since v cc E oc Such dependence has been found in systems where bound-bound transitions 
account for most of the HHG^^''*^. For the case of charge-resonant states of odd-charge molecular ions it has also 
been argued?^ that this simple two-level model can account for the behavior of these molecules in a strong laser 
pulse. A square-root dependence of the cut-off frequency with the intensity of the laser field can be found in several 
experiments^, a fact that is characteristic of a driven two-level system and which might point in the direction of an 
effective two-level dynamics for some systems^. 

HHG and DL are the two mayor features of a driven two-level system with a great range of applications in atomic 
and solid state systems where coherent sources of radiation and control of quantum states are sought. Through this 
work we have established a strong connection between this two phenomena and provided a detailed description and 
a deeper understanding of this system. 
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From the definition of the Floquet-Green operator, we have shown in Ea. (|10|l that the residue of the function 
Qofi{E) at a pole E = e + p is the operator (2x2 matrix) \(t)p){(t)p\, where j^p is the pth Fourier component of the 
Floquet eigenstate corresponding to the quasi-energy e. If we write \(jfp) = Kp\a) + Zp\b), then, in principle, from 
the residue (a matrix) of Go,o we could determine the components Kp, Zp. This is only true, however, provided the 
residue does not turn out to be a diagonal matrix. Unfortunately, this is the case for the system we are considering, 
since, as we mentioned before, Gofi{E) is a diagonal matrix (Veff{E) is diagonal) and therefore, its residue only gives 
us information about \Kp\'^ and \Zp\'^. From Eq.® we can see that the off-diagonal components Gm,o can give us the 
needed components: 



From this we get that at the pole E = (and therefore p = 0), the residue is \4>m) {'Pq'' \- From all these residues for 
all values of m, one could therefore extract all the eigenvector components \(j)m)- 
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Appendix A: components of the FES 




(52) 
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The Gm,o operators can be obtained from Go.o, as it was shown in Ref:~-. The result is 



3=0 

and where (showing exphcitly the energy dependence) 



F+{E) ^ (VT,.o{E))-W = V-XsiE + j) 



for m > 0, 



HiE+j) 

HiE + j) 



and similarly, 



(53) 



(54) 



G-m.o{E) — J-_„^^iG-m+l,Q — ^-jGofi 

j=0 



for m > 0, 



(55) 



with 



FIJE) ^ {VT^,^o{E))-'V = V-XlfiE~j) = 



K(^-J) 

[lv-{E-j) 



It is not difficult to show using Eqs.lO, lf55|l . and that 



and 



(56) 



(57) 



with £^ = ±e. 



As we already mentioned, at E = e the residue of Go,oiE) is proportional to \b) and at E = —e the residue is 
proportional to \a). From this, and using Ea. (|57|l together with Eas. (|54|) and (|56|l . we can write explicitly the two 
Floquet eigenstates of this system: 



p2iLjt iujt 

|r^(t))=7V[....+ -^v^{-e-l)v-{-e)\a) + —v-{-e)\b) + + 



-2iujt 



+ vZ{-e)\b) + -^v+i~e + l)v+i-e)\a) + 



and 



''(t))=N\....^ -^^v-{e-\)vl{e)\h) + —vl{e)\a) + \h) 



+ ^v+{e)\a) + —^v+{e + l)v+{s)\b) + 



(58) 



with N being a normalization constant. If we use the symmetry relations, i^E) = —v'^{E) and Vj^{—E) = —v^ {E) 
we can write this two eigenstates in a more compact way. 



\r'{t))^K{t)\a) + Z{t)\h) 
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and 

\cf>^t))=-Z{tr\a) + K{tr\b) , (59) 

with 

K{t) = N{.... + -^v^{-e - ^)v-{^e - 2)v^{-e - l)v-{-e) + -^v^{-e - l)v-{-e) + 1 
+ + + ^;^<(-e + 3>+(-£ + 2K+(-£ + l)^;+(-£)) 

OO 

= iV[ ^ e-2»-*if2„(-£)] , 

n= — CSO 

Z{t)=N[...—^v-i-e-2)v^ie-l)v-i-e)+ —v-{-s) + -^v+i-e) + —^v+{-e + 2)v+{-e+l)v+{-e)....] 

OO 

= N[Y^ e-^(2"-iMz2„_i(-£)] 

n— — OO 

(60) 

where iV is a normalization constant so that |iir(t)p + |^(t)P = 1, and we have 



K2n{-S) = ■ 
■^2n+l( — e) = • 

and we used sg{x) = x/\x\. 



1 for n = 0, 

ngi ^vf^'H-^ + S9{n){2j - l)yj^''\-s + sg{n){2j - 2))) for \n\ > 

lv+i-£ + 2n)K2n{-s) for2n + l>0, 

iu- (-£ + 2n + 2)K2n+2{-e) for 2n + 1< 0, 
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FIG. 2: Quasi-energy as a function of driving amplitude, for different values of wo > 1 (low frequencies). The functions 
±0.05 Jo (4t;) are plotted with dotted lines for comparison. 
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FIG. 3: Solutions of Eq. (20) as a function of driving amplitude for the resonant case {ujq = 1). The avoided crossings 
structure can be accounted for, in an approximate way, by a set of crossings of the eigenvalues corresponding to the first-order 
FES indicated on the left of the figure. Avoided crossings A and C are mainly due to a 3 photon resonance between the 
corresponding FES; B corresponds to a 5 photon process. 
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FIG. 4: Solutions of truncated versions of Eq. (22): = v^, for different values of n. In a) n=0; b) n=l; c) n=2; d)n=3. 

The labels h — 2n and a — 2m — 1 correspond to terms in the CFs in Eq. (23) of the form 5 + 2n and 5 + 2m + 1 + (Do . As more 
terms are included new avoided crossings appear. The ± signs are used to indicate each one of the two solutions that originate 
at the points 5 = 2n in the y-axis. 
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FIG. 6: Components of the FES for iDq = 1 and for eight different values of the driving amphtude. The logaritmic plots show 
the "plateau" structure of the FES, which gives rise to the well-known plateau studied in HHG. The cut-off location (Ne) 
follows the approximate rule N,; = 2v ~ 1. 
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FIG. 7: The upper panel shows the functions cos^(S(u)) and sin'^{d{v)), for ujo — 0.5. At t = 0, the FES is given by 
10"^ (0)) = cos9{v)\a) + sm6{v)\b). The lower panel shows the functions 9{v) (solid line) and e(«) (dash line). The dash-dot 
Une corresponds to the Struve function d>o^//o(4i;). Points A, B, C, D, correspond to the values of v for which DL occurs. 




FIG. 8: Same as in Fig. 7. Here iuo ~ 1. 
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FIG. 9: Same as in Figs. 7,8. Here ujq = 2. 




FIG. 10: The probability to remain in the initial state |1) as a function of time (r = uit), for 9 different values of the driving 
amplitude and for ujq = 0.5. The lower panels show that the function cos^(2S(w)), at the DL points (A,B,C,D), gives a good 
approximation to the amount of localization Pi in this system. Pi — min(|tiii(r)|^). In the lower left panel the function 
cos^ (2(1)0 J 7^0(41;)) is shown with a dash-dot line for comparison. 
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FIG. 12: As in Figs. 10, 11. Here Cjq = 2. 
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FIG. 13: The emission spectrum of the system, for the "Optical" Initial Condition (OIC) and for two values of the driving 
amplitude: v = 4.5 corresponds to a DL point and v = 4.15 gives ^(4.15) = 0. The Fourier components of are connected 

with a continuous line for better visualization. 
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FIG. 14: Same as in Fig. 13. Here for "Tunneling" initial condition (TIC). 




FIG. 15: Same as in Fig. 13,14. Here for "Complex Tunneling" Initial Condition (CTIC). 



